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Magnetic resonance method for forming a fast dynamic imaging 



The invention relates to a magnetic resonance method for forming a fast 
dynamic imaging form a plurality of signals acquired by an array of multiple sensors 
according to the preamble of claim 1 . The invention further relates to a magnetic resonance 
imaging apparatus for obtaining a fast dynamic image according to the preamble of claim 1 1 
5 and to a computer program product according to the preamble of claim 12. 

In magnetic resonance imaging there is a general tendency to obtain 
acceptable images in shorter time. For this reason there have been developed recently the 
sensitivity encoding method called "SENSE" by the Institute of Biomedical Engineering and 

10 Medical Informations, University and ETH Zurich, Switzerland. The SENSE method is based 
on an algorithm which directly acts on the image as detected by the coils of the magnetic 
resonance apparatus, wherein subsequent encoding steps can be skipped and thus an 
acceleration of the signal acquisition for imaging of a factor two to three can be obtained. 
Crucial for the SENSE method is the knowledge of the sensitivity of the coils, which are 

15 arrange in so called sensitivity maps. In order to accelerate this method there are proposals to 
use raw sensitivity maps which can be obtained through division by either the "sum-of- 
squares" of the single coil references or by an optional body coil reference (see e.g. K. 
Pruessmann et al. in Proc. ISMRM, 1998, abstracts p. 579, 799, 803 and 2087. 

20 The SENSE method is preferred for acceleration of the signal acquisition for 

magnetic resonance imaging, which results in an enormous reduction in operating time. 
However, the method can only be used properly if the coil sensitivity is exactly known. 
Otherwise imperfections will cause fold-over artefacts which provide wrong images. In 
practice the coil sensitivity cannot be estimated perfectly and will be dependent on 

25 fluctuations in time (movement of the patient, temperature influences, etc.). 

Another important problem of the SENSE method is the spatially varying 
noise level in the resulting image. More specifically, the resulting image can have regions of 
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extremely high noise level, due to local "underdetermination" of the information provided by 
the coil patterns. 

It is an object of the present invention to reduce the noise level over the full 
5 image during fast dynamic imaging to a great extent. 

This and other objects of the invention are accomplished by a method as 
defined in claim 1 , by an apparatus as defined in claim 1 1 and by a computer program 
product as defined in claim 12. 

10 

The main issue of the present invention is based on the solution to acquire the 
low-spatial frequencies, i.e. the low-order k-lines, with full density by using the normal 
magnetic resonance imaging, and to acquire the remaining high-spatial frequencies, i.e. the 
high-order k-lines, with reduced density, i.e. with the SENSE imaging method. In more 
15 general terms: it is proposed to have a k-space coverage that is varying over k-space, the 
density of the acquired information being higher in the centre than for the high-spacial 
frequency information. 

An important additional advantage of the imaging in full density of the centre 
20 of k-space is that the SENSE acquisition can be reconstructed without a-priori knowledge of 
coil sensitivity maps. Thus, more accurate images with an improved signal-to-noise ratio can 
be made at the same rate as with the SENSE method. 

These and other advantages of the invention are disclosed in the dependent 
25 claims and in the following description in which an exemplified embodiment of the invention 
is described with respect to the accompanying drawings. It shows 

These and further advantages of the invention are disclosed in the dependent 
claims and in the following description in which an exemplified embodiment of the invention 
30 is described with respect to the accompanying drawings. It shows 

Fig. 1. the regular scanning of k-space according to the well-known 
acquisition scheme, 

Fig. 2, the regular scanning according to SENSE acquisition, and 
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Fig. 3, the acquisition scheme according to the present invention. 

Specific numbers dedicated to elements defined with respect to a particular 
figure will be used consistently in all figures if not mentioned otherwise. 

5 

In figure 1 the regular acquisition scheme according to the well-known normal 
magnetic resonance methods is depicted. In figure 2 the acquisition scheme according to the 
SENSE method is shown. The resulting raw images consist of N mutually overlapping strips. 
Yet, this situation has been measured simultaneously with M coils, whereas M \i N. In fact, 

10 the actual measurement is acquired by using a step in k-space that is N times too large 
compared to the k-space step that Nyquist would prescribe for the regular acquisition of 
figure 1. A typical value of N, also called SENSE factor, is 3, which results in the acquisition 
scheme of figure 2, the data thereof is Fourier transformed. In effect, there exist for each 
pixel of the required image M equations and N unknowns, which can be solved as a set of 

1 5 linear equations. This results in the required image according to the SENSE acquisition 
scheme. 

The main problem of this procedure is that, locally, the set of equations can 
become badly conditioned, i.e. almost singular. Since the input data or coil signals always 
20 contain noise, the ill-conditioning of the set of equations leads to an enormous increase in 

noise in the resulting image. In theory, for a location where the set is completely singular, the 
resulting noise (i.e. the noise in that region of the reconstructed image) becomes infinite. In * 
practice, this is noticed as ugly "noisy patches" in the SENSE images. An additional problem 
is the necessity of coil-sensitivity maps, which are relatively time-consuming. 

25 

In order to speed-up the acquisition with higher accuracy of the reconstructed 
image, it is now proposed to segment the k-space into regions - in a basic set-up into two 
regions. The SENSE factor is different between different regions. Most notably, it is lower 
for the low spatial frequencies. In the most basic set-up the SENSE-factor equals 1 (i.e. 
30 normal acquisition) for the lowest spatial frequencies. This is shown in figure 3. 



From such an acquisition two subsets can be extracted: 
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1) data from a regular SENSE scan, by dropping a number of central lines, in 
order to get the regularity of figure 2, the so-called "S-data". A reconstruction according to 
SENSE would result in a so-called "S-image". 

2) data from a low-spatial frequency scan by dropping all data outside the central 
5 region, the so-called "L-data". This can be reconstructed in a regular manner to a low- 
resolution image, i.e. an image with a low refreshing rate. 

On the condition that the L-data has been acquired with SENSE-factor = 1, 
that data allows the reconstruction of coil sensitivity maps. This means, that a regular SENSE 

10 reconstruction of the S-data can be performed, even if no pre-acquired coil-sensitivity data is 
available. Moreover, the L-data has been acquired with no SENSE or with a moderate 
SENSE-factor, e.g. 1,5. This allows to reconstruct (either by plain FFT or by SENSE) a low- 
resolution image, the so-called "L-image", that can be regarded as stable in terms of signal- 
to-noise ratio. The basic idea is to use the data form the L-image wherever the S-image 

1 5 would be exceedingly noisy. This can be implemented in different ways. 

A first reconstruction set-up, the so-called "hard substitution", reconstructs 
both an S-image and an L-image, and selects the most appropriate one of the two. The 
reconstruction process will "know" where this has to be done: if D is the determinant of the 

20 set of equations to be solved for each pixel of the image. The figure of D strongly varies over 
the whole image, and noise problems are expected in those regions of the S-image where the 
value of D is very low. On the basis of this knowledge the simplest procedure for 
combination of the S-images and the L-images will be, that for every pixel of the resulting 
image the corresponding value of the S-image is selected, if D is above a predetermined 

25 threshold. Otherwise, the corresponding value of the L-image is selected. A slight 

imperfection of this combination technique may be the adverse visual effect of "cut lines" in 
the image, separating low-noise regions from high-resolution regions, i.e. regions with a high 
refreshing rate of the image. 

30 A second reconstruction set-up is the so-called "soft combination", which is a 

variation of the first one: the resulting image I is combined by the following equation: 



I = f(D)L + (l-AD))S 
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where / (D) is a function that approaches 1 for low values of determinant D and approaches 
0 for high D values. L and S are the values of the L-images and the S-images, respectively. 

A third reconstruction method is the so-called "algebraic combination", which 
5 can be described as follows: add to the set of equations, i.e. the basis of the reconstruction of 
the S-image, an additional set of equations, expressing that the resulting image pixel-values 
are somewhat close to the L-values of the L-image. The basis for the algebraic combination 
can be explained in matrix calculus as follows: 

10 The regular SENSE reconstruction is, in essence, solving the following set of 

equations: 

S • p — ffi 

where S is the sensitivity matrix,, so Sy is the sensitivity of coil i on pixel position j. It is an 
M by N matrix, where M is the number of coils, i.e. the number of equations, and N the 
SENSE-factor, or the number of signals folded on each other, i.e. the number of unknowns. 
15 m ("measured value") is a vector of length M: the folded information for every coil; p 
("pixel") is a vector or set of N unfolded pixels. 

The least-square inversion of the problem can be formulated as 
p = (S h • X F~ I ■ S)~ x • S h • • m 

20 where "h" stands for hermitian transposed and *F denotes the M x M noise covariance 
matrix. So far for regular SENSE. 

Now we add an additional set of equations: 
p = fh L 

25 where fh L is the corresponding pixel in the L-image. 

The combined set of equations can be written as 



~ S~ 
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where I N stands for N x N identity matrix. The solution thereof is given by 
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The "extended noise covariance matrix ^P^has to be defined as 
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In this expression R is a diagonal "regularisation matrix" of size N x N. The diagonal 
element R u thereof have to be filled in as the variance of the expected noise in the image, 
times a certain factor (which steers the trade-off between "too noisy" and "too much 
5 resolution loss"). 

The solution of the combined set of equations simplifies to 

p = (S h ■T-S + ir 1 )- 1 (S h rh + R-' -m L ) 

10 This is also a kind of a "weighted addition": if the determinant of S dominates 

over R- 1 (which is, in effect, some constant), then the solution approximates the "regular 
SENSE" reconstruction p = (S h • T" 1 • S)~ x • (S h • • m) ; if the determinant approaches to 
zero, then p « (R~ [ )~ l • (R~ ] • m L ) ~ fh L 

15 The fourth reconstruction scheme is the so-called "full usage". In this scheme 

any of the previous reconstruction methods can be combined with the following procedure: 

reconstruct a combined SENSE-image according to the previous three 
combination methods, 

make an inverse Fourier transformation thereof in the SENSE direction(s), 
20 - substitute the central part of the data thereof with the corresponding form of 

the L-data (note that this "hard substitution" can be replaced by a weighing function which, 
depending on k y , gradually changes from "giving full weight to the L-data" to "giving full 
weight to the combined SENSE data") 

make again a Fourier transformation into the image domain. 

25 

A further option for above mentioned reconstruction methods is filtering of the 
L-data before FFT, by multiplying it with a tapering function, e.g. a Riesz function. In 
addition the L-data can be used for the reconstruction of the coil-sensitivity maps. In regular 
SENSE the coil -sensitivity maps are usually derived by comparison of coil element data to a 
30 body-coil acquisition. However, according to above mentioned methods the body-coil data is 
not or not necessarily available. It is already proposed to take the root-sum-of-squares 
(RSSQ) combination of the coil element data instead. Yet, that RSSQ is devoid of phase 
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information, thus differences with pixels of the L-image may be introduced. Therefore it is 
suggested to take a plain sum of the L-data over the coil elements, which will give a more 
proper result. If the plain sum can be taken after a very basic phase correction, e.g. correcting 
L-data of every coil element to have zero average phase. 

5 

SENSE can also be applied in three-dimensional acquisitions. That can be 
done by a reduced k-line density in k y or k z or both directions, as long as the product of the 
reduction factors does not exceed the number of coils. For the L-data region of k-space, i.e. 
the region of low spatial frequencies acquired with a low SENSE-factor, in principle all 

10 combinations of limitations are possible: limitations thereof in k y or in k z or in both 

directions. The last option is only of interest if there is also a SENSE-reduction in both 
dimensions. Then, the low SENSE-factor region can be e.g. square, circular or rectangular. In 
the case of dynamic scans, is may be that it is not strictly necessary to acquire the L-data on 
every frame: in favourable circumstances that data is only used for rather limited image 

15 areas. Therefore, the L-data of "old" time frames, e.g. of the first dynamic scan, may also be 
sufficient for subsequent images. 



Although the present method is described in respect with SENSE, it can also 
be used for the so-called SMASH acquisition as described in details in Magnetic Resonance 
20 in Medicine, 38, pages 591603 (1997). The main difference with SENSE is that SMASH 
manipulates the image in the Fourier-domain or k-space. 



